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Abstract 

The nonequilibrium total dielectric function lends itself to a simple and gen- 
eral method for calculating the inelastic collision term in the electron Boltz- 
mann equation for scattering from a coupled mode system. Useful applica- 
tions include scattering from plasmon-polar phonon hybrid modes in modula- 
tion doped semiconductor structures. This paper presents numerical methods 
for including inelastic scattering at momentum-dependent hybrid phonon fre- 
quencies in the low-field Boltzmann equation for two-dimensional electrons 
coupled to bulk phonons. Results for mobility in GaAs show that the influ- 
ence of mode coupling and dynamical screening on electron scattering from 
polar optical phonons is stronger for two dimensional electrons than was pre- 
viously found for the three dimensional case. 

I. INTRODUCTION 

The importance of hybridization of collective modes has been pointed out in the context 
of several semiconductor systems. Examples include hybrid optical modes in thin semicon- 
ductor layers coupled intraband-interband excitations of quasi-one-dimensional electron 
systems 0, and, more commonly, plasmon-phonon coupled modes in doped polar semicon- 



ductors |^-|5[ . To the extent that electron scattering depends on the energies and interaction 
strengths of the system modes, a rehable theory of electron transport depends on treating 
mode-coupling effects accurately. 

In the Born approximation, the interaction between a conduction electron and a cou- 
pled electron-phonon system can be represented as an effective electron-electron interaction 
screened by the nonequilibrium total dynamic dielectric function eT{q,uj), which includes 
contributions from both electrons and phonons. Recent work [^|| shows that eT{q,uj) pro- 
vides a systematic way to determine the collision term in the electron Boltzmann equation 
for scattering against dynamically screened coupled electron-phonon modes. The result is 
the sum of an electron-electron collision term and an electron-LO phonon collision term 
that includes plasmon-phonon mode coupling. Both interactions are dynamically screened 
by only the electronic part of the total dielectric function for the electron-phonon system. In 
the random-phase approximation (RPA), the electron-phonon part contains a phonon self- 
energy that arises from the polarization of the electron gas. The self-energy correction 
modifies the longitudinal-optical (LO)-phonon dispersion in doped polar semiconductors, 
producing hybrid normal modes with phonon strength in each. Numerical results |^ for 
bulk n-type GaAs show that mode-coupling and dynamical screening should significantly 
influence electron mobilities in modulation doped structures. 

The present paper describes numerical methods for exactly solving the low-field Boltz- 
mann equation for two-dimensional electrons coupled to bulk LO phonons, including dy- 
namical screening and mode coupling for arbitrary electron degeneracy and spherical energy 
surfaces. In particular, the method for including phonon dispersion in the collision integrals 
is discussed. The phonon distribution is approximated by its equilibrium form, and the 
plasmon-pole approximation is used in the LO phonon self-energy to determine the hy- 
brid mode frequencies. Results for GaAs show that mode coupling and dynamical screening 
are more important for two dimensional electrons than for the three dimensional case. 
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II. INELASTIC SCATTERING FROM HYBRID PHONONS WITH DISPERSION 



The collision integrals in the electron Boltzmann equation for doped polar semiconductors 
contain the momentum- dependent frequencies of the hybrid LO phonon modes 0. Using 
the plasmon-pole approximation in the phonon self-energy, the differential scattering rate 
due to transfer of momentum q = k — p to the hybrid phonons is 

— 2M^ 

Ty^O(k,p) = -— ^-—[N{uk,p) + l]lm[D+{q,u,,,) + ^-(g,^,,,)] (1) 
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where M| = Vg{e^ — eQ^)hujlQ/2ujTo and hcukp = Ek — Ep. Here, ojlo and ujto are the 
LO and transverse optical phonon frequencies, while eoo and eo are the high-frequency and 
static dielectric constants, respectively. In two dimensions, Vq = 27re^/geoo- The frequencies 
Up, and uj_ are given by 
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where ujp = {2Tme'^ q / m* eooY ^'^ is the plasmon frequency of the two-dimensional electron gas 
with concentration n and effective mass m*. The weight factors {uj\ — uj'^){uj']^ — uj'^_)~^ in 
Im[D^] give the phonon strength in each of the hybrid uj± modes, so that the differential 
scattering rate W^^ is the rate for scattering from only the phonon component of the hybrid 
modes. The screening function in (|l]) is the temperature-dependent RPA dielectric function, 
e{q,uj) = 1 — v'^P{q,u!), for the two-dimensional electron gas. 

The iterative procedure used here for solving the low-field Boltzmann equation was de- 
scribed in reference [^. One assumes a linear form for the nonequilibrium electron distri- 
bution, /(k) = fk + XkQk, where /° is the equilibrium (Fermi-Dirac) distribution, Xk is the 
cosine of the angle between the electric field F and k, and is an unknown function that 
is linear in F. Neglecting all scattering mechanisms except electron-LO-phonon gives 
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^Lo[9] = -E^?P^'cp{l^^°(p,k)[l-/°] + iy^O(k,p)/°}. 
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and Xkp is the cosine of the angle between k and p. 

The form of the phonon Green's function in equation (2) imphes that the energy conser- 
vation relation between initial and final electron states, Ek — -Ek-q = i^t^lo'), is dependent 
on the momentum transfer q through the phonon dispersion. In three dimensions it is 
convenient to choose q as an integration variable when calculating z/^o aiid t^q- In this 
case, one must determine the integration limits by finding the points where the phonon 
dispersion curve oj{q) intersects the absorption region q^ + 2kq < 2m*uj{q) /% < q^ — Ikq 
and the emission region 2m*uj{q) < —q^ + 2kq. For the case of two-dimensional electrons, 
it appears that q is not a preferred choice for integration variable, since the integrand then 
diverges at the integration limits. Instead one can choose d^p, the angle between k and p. 
In this case, one must determine the implicit relation between q and 6kp at each integration 
step by finding the (nonzero) root of the function /(g) = — A;^ — + 2kp± cos 9kp with 
p± = ±2m*uj±/hy/\ 

Scattering from the LO phonon part of the hybrid modes gives the inverse lifetime TEo{k) 
and functional i'loWi^ 

-rEoik) = E a| r de,pI{q,uJx)\N{uJx) + 1 - fJ&iEk - Tiu^) 
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where k^ = (/c^ ± 2m*ux/fiY^'^ and 
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III. RESULTS AND CONCLUSIONS 



Figures 1 and 2 present electron mobility as a function of concentration in GaAs at 300 
K and 77 K, showing the effects of dynamic screening and mode coupling on the electron-LO 
phonon interaction. The mobilities were calculated by averaging the electron velocity over 
the nonequilibrium distribution /(k), found by solving the iterative equation (5). Without 
plasmon-phonon coupling, dynamic RPA screening gives mobilities that are lower than static 
RPA or Thomas-Fermi results, but higher than the unscreened case. The difference between 
static and dynamic screening has a strong momentum dependence. The inverse dynamic 
RPA dielectric function evaluated at = ujlq has a strong peak (e^^ > 1) at small q where 
the plasmon dispersion crosses ujlO) but is smaller than the inverse static dielectric function 
at larger q values. The small q values, where antiscreening and enhanced forward scattering 
occur, are accessible only for high energy electrons, since the minimum q value allowed by 
energy conservation is qmin = \k — [k"^ + 2m*ujLo/K)^^'^\- Previous authors have noted 
similar effects in calculations of dynamical screening of the electron-LO phonon interaction. 

When plasmon-phonon coupling is included in the plasmon-pole model, the phonon spec- 
tral function — 7r^"'^Im[D(g, cj)] has two peaks with frequencies uj+ and uj^ given by (^. 
Mobilities are expected to be lower when mode coupling is included, because of increased 
low-energy electron-phonon scattering due to the uo_ mode. Figures 1 and 2 show that this 
is in fact the case. As in three dimensions P|, mode coupling works to reduce electron 
mobility especially for low densities at 77K where the thermal occupation of the low-energy 
hybrid mode is exponentially larger than the high-energy hybrid mode or uncoupled phonon 
mode. The effect is more pronounced in two dimensions, consistent with previous conclusions 
concerning the effect of dimensionality on polaronic damping. 
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FIGURES 

FIG. 1. Effects of screening and mode coupling on LO phonon limited mobility of 

two-dimensional electrons in GaAs at 300K. The dynamically screened coupled mode mobility calu- 
lation (solid curve) is compared to mobilities determined by scattering from uncoupled LO phonons 
screened in the Thomas-Fermi approximation (long-dashed curve), static RPA (short-dashed curve), 
dynamic RPA (dotted curve) , and unscreened (dot-dashed curve) . All other scattering mechanisms 
are neglected. 

FIG. 2. Effects of screening and mode coupling on LO phonon limited mobility of 
two-dimensional electrons in GaAs at 77K. Mobility curves are represented as in Figure L 
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